Phase diagrams of XX Z model on a depleted square lattice 
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Using quantum Monte Carlo (QMC) simulations and a mean field theory, we investigate the 
spin- 1/2 XX Z model with nearest-neighbor interactions on a periodic depleted square lattice. In 
particular, we present results for 1/4 depleted lattice in an applied magnetic field and investigate 
the effect of depletion on the ground state. The ground-state phase diagram is found to include 
an antiferromagnetic (AF) phase of magnetization m z = ±1/6 and an in- plane ferromagnetic (FM) 
phase with finite spin stiffness. The agreement between the QMC simulations and the mean field 
theory based on resonating trimers suggests the AF phase and in-plane FM phase can be interpreted 
as a Mott insulator and superfluid of trimer states, respectively. While the thermal transitions 
of the in-plane FM phase are well described by the Kosterlitz-Thouless transition, the quantum- 
phase transition from the AF phase to in-plane FM phase undergo a direct second-order insulator- 
superfluid transition upon increasing magnetic field. 

PACS numbers: 75.10.Jm, 75.45.+J, 05.30.Jp, 75.40.Mg 



I. INTRODUCTION 

It is well known that depletion would weaken the 
long-range spin correlation in various quantum-spin 
systems^ - — Reduced dimension arising from depletion 
would enhance quantum fluctuations in the system that 
consequently diminishes the spin-spin correlation in the 
long-range regime. It is known that the zero-point fluctu- 
ations in one-dimensional spin-1/2 Heisenberg model lead 
to power law decay in the spin-spin correlation and there- 
fore destabilize the classic Neel state. Depletion of sites in 
two-dimensional magnetic systems would, therefore, ulti- 
mately destroy the long-range order (LRO) in the ground 
state by lowering dimensionality. One example is the de- 
struction of the LRO of triangular spin-1/2 Heisenberg 
antiferromagnet by depleting 1 /4 of the spins that leads 
to a Kagome latticed Similar phenomenon also occurs in 
triangular Bose-Hubbard model. The LRO phase in the 
triangular lattice^ appears when nearest-neighbor (nn) 
repulsive coupling V > 5.13 but can only be stabilized 
for larger V > 8 in the Kagome lattice^ because of the 
enhanced quantum fluctuations caused by the site de- 
pletion. The LRO phase also emerges as an interesting 
valence-bond solid where bosons are delocalized around a 
subset of hexagons in the depleted system. Unlike these 
frustrated systems, another experimental important ex- 
ample is the 1/5 depleted Heisenberg model of CaV^g 
where the increased quantum fluctuations are not suffi- 
cient to destroy the LRO^^ Heisenberg model of 1/4 
depleted square lattice has also been investigated^ Al- 
though spin correlations are weakened due to the deple- 
tion, quantum fluctuations, again, are not strong enough 
to destabilize the LRO. Given its simple lattice structure, 
it is surprising to note that a complete investigation of 
the XX Z model in the 1/4 depleted square lattice is still 
missing. 

On the other hand, intense efforts have recently been 
focused on the exotic supersolid 10,11 phase on a variety 
of models£ii2~— . One natural question to ask is whether 
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FIG. 1: The unit cell and site indices of the 1/4 depleted 
square lattice. The spin pattern represents the m z = 1/6 AF 
phase. 



depletion would, via reducing diagonal LRO, allow the 
possibility of coexistence of the off-diagonal and diagonal 
LROs. Therefore, it will be worthy to address this issue 
in the simplest depleted system. It is the object of this 
work to investigate the role of depletion on the ground- 
state phase diagram of the 1/4 depleted XX Z model 
by employing both the quantum Monte Carlo (QMC) 
simulation and the timer basis mean field (MF) theory. 



II. MODEL 

The lattice structure and the unit cell of the system are 
shown in Fig[TJ Each unit cell contains three spins with 
two of which (sites 1 and 3) are equivalent by symmetry. 
As shown in the figure, site 2 connects to sites 1 and 3 and 
has the coordination number z = 4, while for the sites 
1 and 3 z = 2. This difference in coordination number 
z leads to distinct chemical potentials fi on these sites 
when transformed the XX Z spin model into the hardcore 
Bose-Hubbard model. 

In this study, we consider the standard spin-1/2 XX Z 



Hamiltonian: 
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AS*S? + -(S+Sr + h.c.) - /i£ S, 



(1) 

in the depleted square lattice. Here the common nota- 
tions of spin operators S z , 5 ± are used. To facilitate fur- 
ther analysis, Hxxz c & n be decomposed into the inter- 
cell and intra-cell parts and is written as: 



H s — ^ (H ai + Hb, 



(2) 



A^(^ + S? 3 ) + -S+(S a + S~ 3 ) + h.c. 
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+/i.c. 



where 7J a i and are the Hamiltonians responsible for 
the inter-cell and intra-cell couplings, respectively. 
now sums over all unit cells and YJ Q sums over the sites 
within a cell. H s is equivalent to a hardcore boson Hub- 
bard model under the transformations SZ, — bf n and 
SL = n, 



s. such that: 



W^+i^ + ~(A + 3/i) 



(3) 



Vni 2 ("ii + ^3) + -^>f 2 {b iX + b i3 ) + h.c. 



Vn l2 (n l+£1 + n i+ y 3 ) + -bj 2 {b l+il + b l+ ~ 3 ) 
+h.c, 



in which fi\ = fi 3 = h + A and [i% = h + 2A due to 
the distinct coordination number 2 for the two types of 
sites. Hereafter, we will mainly focus on the discussion 
of the spin model but will also use the boson language if 
appropriate. 



III. QUANTUM MONTE CARLO 

The spin XX Z Hamiltonian can efficiently be simu- 
lated by the standard stochastic series expansion (SSE) 
(Ref. 17) approach without difficulties. Most of the re- 
sults are done on lattices of 24x24 unit cells and up to 
36x36 unit cells are considered in the study of quantum- 
phase transitions. Temperatures are chosen to be in- 
versely proportional to the lattice size such that T oc 
X/(2L). The essential order parameters include the spin 
stiffness p s (or superfluidity in boson language) that sig- 
nals the in-plane FM phase. It is related to the winding 
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FIG. 2: (Color online) The ground-state phase diagram ob- 
tained by the QMC (denoted by black dots) simulations of 
lattice size 28x28. AF phase, characterized by the magneti- 
zation plateau m z = ±1/6, extends further into the in-plane 
FM phase as explained in the text. Results of a trimer-basis 
MF theory (blue dashed line) agree very well to the QMC 
data. For A < 1, phase separation (PS) is observed at h — 
while the ground state is a uniform superfluid for A > 1. FP 
denotes the fully polarized states for large fields. The vertical 
red line shows the scan for the order parameters in the Fig[3] 



number fluctuation of the updated loops in the SSE algo- 
rithm and can be easily implemented.— The broken sym- 
metry of the lattice structure is reflected by the measured 
spin-structure factor 



S(Q)/N 
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in the simulation. We set the lattice constant of each 
unit cell as unity and focus only on the wave vectors 
Q = (2?r, 2tt), (2tt, 0) and (0, 2vr) by which the lattice 
symmetry of the depleted lattice is characterized. 



IV. GROUND-STATE PHASE DIAGRAM 

The ground-state phase diagram of hardcore Bose- 
Hubbard model in the undepleted square lattice^ con- 
tains simply a lobe of half-filled solid of checkerboard 
structure (m z — Neel state) surrounded by a super- 
fluid phase. Upon depletion of 1/4 sites, the m z = 
Neel phase splits into two antiferromagnetic (AF) phases 
with magnetization m z = ±1/6 due to one missing spin 
in each unit cell as illustrated in the Fig. [1] As a con- 
sequence, the original m z — lobe in the ground-state 
phase diagram now splits into two AF lobes in the de- 
pletion system as shown in the Fig. [2] Due to the ge- 
ometrical frustration at this magnetization, there is no 
m z = state at zero field, and instead, a phase separa- 
tion is observed for A > 1. For smaller A, however, the 
zero-field ground state is a uniform in-plane FM charac- 
terized by a finite spin stiffness p s . In contrast to the 
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fleet the depletion-induced internal structure of the unit 
cells. The superfluid indeed does not break the crystal 
symmetry. This differs from the usual spin supersolid 
where the crystal symmetry is broken in addition to the 
broken gauge symmetry. Therefore the regular depletion 
in our case has not lead to any supersolid phase in the 
ground-state phase diagram. 

At zero field, all sublattice magnetizations are the same 
and it indicates a uniform FM or superfluid phase. To 
shed light on the ground-state wave function, we con- 
struct a simple trimer-basis MF theory which surpris- 
ingly, agrees very well to the numerical results. 



V. MEAN-FIELD THEORY 



FIG. 3: (Color online) QMC data of order parameters of the 
(a) total and sublattice magnetizations, (b) spin stiffness p a , 
and (c) structure factors as a function of field h for A = 0.8. 
The dashed lines in (a) denote the corresponding mean-field 
magnetizations. 

undepleted lattice, in which the AF LRO occurs only for 
A ^ I, the AF phase in Fig. [2] extends further into the 
A < 1 regime and leads to the re-entrance of the in-plane 
FM phase when scanning along the increasing h (see Fig. 
[3]). This remarkable result indicates the regular deple- 
tion in our model in effect surpasses the quantum fluctu- 
ations and enhanced the AF LRO for smaller A > 0.57. 
Furthermore, the quantum-phase transition from AF to 
in-plane FM is continuous, in great contrast to the first- 
order phase transition from Neel to FM phases appeared 
in the undepleted square lattice. We stress that these 
observations are different from the case of the related 
Kagome lattice, which is essentially a regularly depleted 
triangular lattice. The region of the solid phase shrinks 
in that case and the solid to superfluid phase transition 
remains first order. Detail discussion of phase transitions 
will be offered in Sec. IVII In FigEJwe present the order 
parameters as a function of h with a typical A = 0.8. 

The average magnetization m z increases monotonically 
as h except within the AF phase which exhibits as a 
plateau of fixed m z — 1/6 as plotted in the Fig. [3] 
(a). To gain further insight, the sublattice magnetiza- 
tions m Q is also measured. The plateau phase has a 
finite-structure factor at Q = (2tt, 2w), which represents 
the checkbroad crystal structure, while the spin stiffness 
p s vanishes. Sublattice magnetizations mi and 777,3 are 
identical but are different from the m.2 due to the dis- 
tinct geometry. Away from the AF plateau phase, p s are 
finite for all h smaller than the critical field h C 3, beyond 
which all spins are fully polarized. However, these in- 
plane FM (superfluid) states are not uniform in space as 
illustrated by the finite structure factors 5(Q ) and the 
distinct sublattice magnetizations. The superfluid con- 
densates vary at site 2 and sites 1 and 3 such that the 
ground state is a non-uniform superfluid. Nevertheless, 
one has to note that the finite values of 5(Q ) only re- 



We begin with the exact ground-state wave function 
of the spin trimers in a unit cell, i.e., ignoring the inter- 
cell coupling Hbi in the original Hamiltonian H s in the 
Eq. [2] Denoting the spin-trimer basis as IS1S2S3) with Sj 
represents the spin at the site i, the four lowest-energy 
eigenstates of the trimers are: 

ip +i = _L_(|tt;) + utt)+«int)) 

2 v 1 + or 
<P+i = I tit) 

a = ~ (A + VA 2 + 8) 
with eigenvalues 
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respectively, and the states i/3_i and with spin up 
and down interchanged. For < h < (A — a)/2, the res- 
onating trimer states ip + 1 and are the lowest energy 
and the first excited state, respectively. Up to now the 
analysis is exact. These states are no longer eigenstates, 
however, once the inter-cell coupling is taken into 
account. To proceed we neglect the higher-energy states 
and construct a superfluid wave function of the trimer 
states ip±i such that: 

*o = I] ( u °¥>+ii + voe*V$i) (5) 

i 

for small external field h. Here uq and vq are related 
by Uq + Vq = 1 and are determined by the minimization 
of the total energy E = (*o|-ffs|*o)< The global phase 
9 will not alter the total energy and is set to zero for 
convenience. If </2-|i i s considered as the vacuum state, 
the creation of a V+ij state is equivalent to adding a 
spin 1 hardcore boson and ^0 represents a superfluid of 
those hardcore bosons. At zero magnetic field h — 0, the 



occupation of both ip±ii states are identical with uq = 
vq = that leads to the zero total magnetization (m z — 
0). Increasing h will raise the number of ip + i i states by 
flipping spins in the ip_ ij states. The condensate density, 
cc uqVq, reduces as uq increases and vq decreases from 
l/v2- Ultimately, «o drops to zero at a critical field 
h c i as represents a Mott insulator of <p+if with a 
spin gap that manifests as the magnetization plateau of 
the AF state in FigO] Higher external fields will lower 
the energy of V+'^i an d finally close up the spin gap at 
h = h C 2- The ground state is then represented by: 



*r=n 



Ul<p, 



+ vie *ip. 



(6) 



as another superfluid (SF) of hardcore bosons ip + s. i - The 
occupation of m which spins are all polarized, 

increases as h, until the system is fully polarized at 
/i C 3 = A — a. The critical fields h c \, h C 2, and h c ^ de- 
fine the MF phase boundaries which are plotted in the 
ground-state phase diagram (Fig. for comparison with 
QMC results. The agreement is good enough to suggest 
that the simple trimer basis MF theory is capable of cap- 
turing the essential low energy physics of the depleted 
system. It is noted that mean field calculations on the 
square and triangular lattice s 19 i 20 also yield similar good 
agreements to the exact numerical results. Nevertheless, 
significant derivations are observed in the vicinity of the 
tip of the AF phase. This is attributed to the proximity 
of the energy levels of the <p_i, <p + i, and if + i i so that 
all of them should be taken into account in the wave 
function in order to achieve a better agreement with the 
QMC results. 

This trimer-basis MF theory predicts all the transitions 
are continuous and of the SF-insulator type, in contrast 
to the Kagome lattice in which discontinuous transitions 
are found. The QMC results of the quantum-phase tran- 
sitions will be addressed in the next section. 



VI. THERMAL AND QUANTUM-PHASE 
TRANSITIONS 

For the case of undepleted square lattice, 12 the 
quantum-phase transitions from the in-plane FM phase 
to the AF phase (superfluid to solid phase in boson lan- 
guage) is apparently discontinuous because of the differ- 
ent broken symmetries in both phases. Similar situation 
occurs in the triangular lattice, 5 albeit a supersolid phase 
exists in some small regimes of the phase diagram such 
that continuous quantum-phase transitions are observed 
from the supersolid phase to the superfluid phase. Away 
from the supersolid phase, direct transitions from the 
superfluid phase to solid phase are of first order. Upon 
depiction, since the superfluid and solid phases break dis- 
tinct symmetries, a weakly first-order phase transition is 
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FIG. 4: (Color online) Finite-size scaling of p s at h c i, the 
first critical field of in-plane FM to AF phase transition in 
Fig. [3] for A = 0.8. (a) Data of p s L 2 for different system 
sizes intersect at a critical field h c i = 0.2568(1). (b) Data of 
p s L 2 collapse well into a universal scaling curve when plotting 
against the (h — h c i)L 1 ' v . The dynamical critical exponent is 
set to z = while correlation exponent v — 0. The tempera- 
ture used is P = L z /12. The same for Fig. [5]and[6] 



recently observed^. These results are in contrast to our 
findings of the depleted model H s . 

To investigate the quantum-phase transitions in detail, 
we performed a finite size scaling (FSS) analysis for the 
continuous AF to in-plane FM phase transitions. In two- 
dimensions, close to a critical point K c , p s scales as2i, 



p. = L-F p .[L 1 l v {K-K ),p/L'\ 



(7) 



where L is the linear system size, v the correlation expo- 
nent, z the dynamical critical exponent, and K c —K is the 
distance to the critical point. Given the scaling relation, 
data of different system sizes L should cross at the tran- 
sition point K c when p s L z is plotted against K as long 
as (5 changes as L z so that the second argument in F Ps 
remains constant. Equation [7] also implies the collapse 
of data onto a universal curve of the scaling function F Pa 
once p s L z is plotted against (K — K^L 1 ^ for the appro- 
priate critical exponent v. For A = 0.8, we performed 
the finite-size scaling of p s at two critical magnetic fields 
h c i and h C 2- As illustrated in Figs. Ufa) and[S](a), data of 
different system sizes intersect, respectively, at the same 
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FIG. 5: (Color online) Finite-size scaling p s at /i C 2, the second 
critical field of AF to in-plane FM phase transition in Fig. [3] 
for A = 0.8. Data also appear to collapse for z — 2 and 
v = 0.5, though not as well as in Fig. [4] for h c i. 



critical points of h cl = 0.2568(1) and h c2 = 1.1426(2) 
for z = 2. Furthermore, data of p s L 2 in Figs. 0Jb) an d 
E{b) appear to collapse well for v = 0.5 as expected from 
the superfluid-insulator universality class^i supports the 
notion of continuous phase transitions between AF and 
in-plane FM phases. 

Similar analysis applies to the melting of AF phase to 
FM phase crosses the critical A c (horizontal line in Fig. 
[SJ. In Fig. [5J the QMC data scale well as stated in Eq. 
[7] for the critical exponents z = 1 and v = 0.67 with the 
critical A c = 0.5685(5). All these results suggest the ob- 
served quantum-phase transitions of AF to in-plane FM 
phases obey the Landau-Ginzburg- Wilson theory in 2D. 
In addition, the observation of a single peak in the his- 
togram of QMC data at the critical points implies the 
possibility of weakly first-order phase transition is un- 
likely. 

The critical behaviors of the depleted lattice can be 
understood as one noticed that there is in fact no broken- 
crystal symmetry in the AF phase. As mentioned in the 
Sec. IV, the finite-structure factors 5(Qo) only reflect 
the internal spin structure within unit cells which origi- 
nated from the two different types of sites in the lattice. 
Although 5(Qo) are found finite in both phases, the crys- 
tal symmetry is indeed preserved. Therefore the phase 
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FIG. 6: (Color online) Finite-size scaling of p 3 crossing the 
tip of the lobe in Fig. [2] (red horizontal line in Fig. [2]) with 
m z = 1/6. Using z = 1 and v = 0.67, (a) the critical A c is 
found to be A c — 0.5685(5) and (b) data collapse well onto 
the universal scaling curve. 



transition is expected to be continuous between the in- 
sulating AF and the in-plane FM phase. On the other 
hand , the situation is different in the Kagome lattice, as 
well as the undepleted square lattice, where all sites are 
equivalent. The observed finite-structure factors indicate 
that the spin structure of the AF phase breaks the crystal 
symmetry of the system. Consequently, the phase tran- 
sition to the in-plane FM phase is expected to be first 
order. This explains the distinct critical behaviors upon 
depletion in these closely related systems. 

Besides the quantum-phase transitions, we also investi- 
gate the thermal phase transitions which are found to be 
well described by the Kosterlitz-Thouless (KT) theory^ 
in 2D. For T > 0, only quasi-long-range order is possible 
in the xy plane as continuous gauge symmetry cannot be 
broken in 2D. Nevertheless the spin stiffness p s can re- 
main finite until reaching the KT temperature T K t, which 
can numerically be determined by noting that the univer- 
sal jump of p s at T KT is given by 



Ps{Twy) — — 2kt- 

7T 



(8) 



6 




FIG. 7: (Color online) p s is plotted as a function of T for 
different system sizes L near the KT transition temperature 
T KT for A = 0.8 at h = 0. The dashed line p s (T) = §T in- 
tersects p 3 curves of different L at temperatures T*(L), which 
follows the logarithmic behavior given by Eq. [9] as shown in 
the inset. 



h=1.8 




FIG. 8: (Color online) Same as Fig. [7]but now h = 1., 



Taking into account the logarithmic finite size correc- 
tions, the intersections of the line p s (T) — —T and the 
QMC data for different system sizes occur at 



work. The phase diagram contains an in-plane FM phase 
and a AF phase, where the AF phase extends further 
than the case in the undepleted lattice due to the regu- 
larity of the depletion. In the in-plane FM phase, in spite 
of the coexistence of finite-structure factor 5(Q ) and 
p s , the ground state does not break the crystal symme- 
try and is not a supersolid state. Quantum-phase tran- 
sitions from the AF to in-plane FM are demonstrated 
to be continuous, in stark contrast to the case of unde- 
pleted lattice. The critical exponents v and z retrieved 
from the FSS are consistent with the prediction from 
the superfluid-insulator universality class. On the other 
hand, numerical data of the thermal phase transitions be- 
long to the KT type as expected for the 2D superfluids. 
Apart from the QMC results, a trimer-basis MF theory 
is proposed to gain more insights of the model. Based 
on the trimer states of each isolated unit cells, the MF 
theory projects out higher-energy states and attempts to 
capture essential physics of the system by retaining the 
lower-energy states. The AF is represented as Yii 
and the FM can be expressed as superfluids of <p + ! in 
the background of or (p + s.. The agreement between 
QMC and MF results is remarkable, given the simplicity 
of the trimer-basis MF theory. 

While a triple laser beam design in cold atom gases2£ 
has been proposed to simulate a Kagome lattice, whether 
similar techniques may apply to the 1/4 depleted square 
lattice is worthy to pursue. Furthermore, this model also 
related to the Bose-Hubbard model with both nn and 
next-nearest-neighbor couplings. Recent study on this 
model found a commensurate SS of density 1/4, which 
has similar structure of the 1/4 depleted lattice studied 
here2^. One interesting issue is if the nun coupling (i.e. 
diagonal bonds in our case) would lead to unexpected 
types of quantum phases. Moreover, adding only half of 
the diagonal bonds in the 1/4 depleted lattice leads to 
the Kagome lattice in which very different ground states, 
including VBS solid, is observed. The investigation of 
how do the two distinct phase diagrams merge as the 
diagonal coupling is increased will also be important for 
the understanding of phase transitions of VBS solids. 



T*(L)=T KT 1 



21n(L/L ) 



(9) 



As shown in Figs. [7]and[5]for h = and 1.8, respectively 
for L up to 36, the intersections T* follows well the log- 
arithmic correction of Eq. [5J For infinite system, p s 
will vanish at the measured T K j and the system becomes 
paramagnetic above Tkt. 



VII. SUMMARY 

A comprehensive study in the ground-state phase di- 
agram of the 1/4 depleted square lattice is given in this 
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